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CN ■ Abstract: In this paper we investigate the relation between the bulk and boundary 

in AdS/CFT. We first discuss the relation between the Poincare and the global 
vacua, and then study various probes of the bulk from the boundary theory point 
of view. We derive expressions for retarded propagators and note that objects in 
free fall look like expanding bubbles in the boundary theory. We also study several 
Yang-Mills theory examples where we investigate thermal screening and confinement 
using propagators. In the case of confinement we also calculate the profile of a flux 
tube and provide an alternative derivation of the tension. 
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1. Introduction 



After the discovery of the AdS/CFT correspondence ]]]] (see also [0, |3j) and its more 
elaborated formulation in [fflH], numerous supergravity calculations have been used 
to learn more about supersymmetric Yang-Mills theories. Particularly interesting is 
the calculation of Wilson loops in || |7| which have also been generalized to finite 
temperature || |9|, [U], [11], [12| . Ultimately, one hopes to be able to use the corre- 
spondence to study also non supersymmetric Yang-Mills theories, preliminary steps 
in this direction have been taken in [13, 14 1. In particular there exist calculations of 



glueball masses that agree reasonably well with lattice calculations ||T5|, [L6], [TTJ, [IS 
It is however far from clear how to connect these models with real QCD in a rigorous 
way [19 



In pi| , p2| the relation between scale in the boundary and radial position in the 
bulk was studied in some detail. This is especially interesting from the point of view 
of holography, i.e. the capability of the boundary theory to fully encode the bulk 
theory. (Note that this can be expected only if the bulk theory includes gravity, while 
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the use of a bulk theory to describe a boundary theory can be expected to be more 
general). In this paper we will continue the investigation of this correspondence. 

In the first section we clarify some issues concerning vacua defined using various 
coordinate systems. Next we derive exact retarded propagators that can be used to 
calculate boundary fields for time dependent bulk configurations. In section three we 
proceed by considering various Yang-Mills theory examples. We derive the screening 
length and the F^ u expectation value at finite temperature using WKB techniques. 
For the confining case we calculate the shape of a flux tube and verify that the 
integrated energy density gives a tension that reproduces the confining potential 
calculated using other means. 



2. The Global and the Poincare Vacuum 

A basic starting point in field theory is the vacuum state. In curved space such as the 
AdS space the appropriate choice is often non-trivial and has a great impact on the 
physical picture. For example, in investigating the thermodynamics of black holes it 
is crucial to impose a correct initial vacuum state to see the Hawking radiation. An 
often used vacuum choice in the AdS/CFT literature is the vacuum corresponding 
to the Poincare coordinate system, the Poincare vacuum. For example, the Poincare 
vacuum is the chosen initial vacuum state in studies of thermodynamics (see e.g. 



23, 24, 25, Eg) of BTZ black holes [Eg, 28]. A slightly puzzling issue is that the 



Poincare coordinates cover only a part of the AdS manifold so there is a coordinate 
horizon. It would seem preferable to use coordinates that cover the whole manifold, 
the global coordinates, and use the associated global vacuum state as a starting 
point. It is therefore important to understand if there is a difference in the two 
vacuum choices, and what the difference will be. 

In this section we study the relation between the global and the Poincare vacuum. 
First we define the vacua in the framework of quantized free bulk fields and the 
corresponding operators in the boundary theory, and make a comparison between the 
states in the Poincare and global coordinates. We then move into a more general level 
and study the Virasoro generators in the Poincare and global coordinates, and relate 
the two vacua into each other using the transformation properties of the Virasoro 
generators. By the bulk/boundary correspondence, it should not matter whether the 
states are compared using their bulk theory or their boundary theory representation, 
since both represent the same Hilbert space. We will therefore mostly focus on the 
representation of the states in the boundary theory. 

2.1 AdS coordinates 

In this section we collect some results concerning various coordinate systems for AdS. 
We will focus on AdS^, but the extension of the results to the general case should be 
obvious. 
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AdS^ can be denned as the submanifold of R 2,2 denned by the equation 

U 2 + V 2 -X 2 -Y 2 = 1 (2.1) 

where U,V,X and Y parameterize the embedding R 2,2 . 

In the following we will use three coordinate systems, namely global (r, /i, #), 
Poincare (t,xo,x) and BTZ (t,<p,r) defined through 

U = cosh fi cos t = - — (1 + x 2 + Xq — t 2 ) = r cosh(r + — r_t) 
2xq 

t 



V = cosh fj, sin r = — = \J r 2 — 1 sinh(— r_0 + r+t) 
x 

X 

X = sinh a sin 9 = — = r sinh(r + — r_t) 
x 

1 



y = sinh/xcos^ = (1 - x 2 - xl + t 2 ) = Vr 2 - 1 cosh(-r_0 + r + t) (2.2) 

2x 

The BTZ coordinates parameterize only a patch within the Poincare space which is 
also only part of the whole AdS^ space spanned by the global coordinates. If the 
identification <fi = (j> + 2tt is made in BTZ coordinates then the space becomes a 
black hole and the parameters r + and r_ represent the radii of the outer and inner 
horizons. 

The boundary of Ad S space is at \i — > oo, xo — > 0, r — >• oo in the different 
coordinates. Near the boundary the coordinate transformations look like conformal 
transformations between boundary variables. It is easy to see that the conformal 
transformations are 

z = tan (^2^ = tanh ^(r + — r -)~7^ 

z = tan = tanh Ur+ + r-)^\ (2.3) 

with 

Z = t + X,W = T + 6,(f) + = t + (f) 

z — t — X , W — T — 9 , (p- — t — (j) 

2.2 Quantization of the Bulk Fields and Boundary Operators 

We begin with a short review of the operator quantization scheme in the bulk and 



boundary theory in Minkowski signature, following the references [2l||[22|. A free 
bulk field is promoted to a quantized operator through an expansion 

$ = Y, a kM r ,x) + h.c. (2.5) 

k 

where <ftk are the normalizable modes in the bulk, in the chosen coordinate system, 
and 0^ are the associated annihilation operators. Here r denotes the radial coordi- 
nate, and x labels the boundary coordinates. 
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The vacuum |0) is denned as the state that satisfies 



flfc|0) = (2.6) 

for all annihilation operators. If we use the Poincare coordinates, we call the result- 
ing vacuum the Poincare vacuum, similarly the global coordinates define the global 
vacuum. 

The boundary operator O coupling to the bulk field $ has a similar operator 
expansion 

0=J2hMx) + h.c. (2.7) 

k 

where (pk are the boundary data extracted from the asymptotic behavior of the 
normalizable bulk modes near the boundary: 

<t>k ~ r 2h +M%) • (2-8) 

By the bulk/boundary correspondence, the vacuum defined by 

b k \0)=0 (2.9) 

corresponds to the same state in the abstract Hilbert space of the theory as the 
vacuum defined using the bulk annihilation operators a&. We have only changed the 
representation of the Hilbert space from bulk to boundary theory language. 

Let us now take a closer look into the boundary theory in the two coordinate 
systems. The boundary is the boundary of the covering space of AdS 3 (in this section 
we will focus on 2+1 dimensions). It is an infinite cylinder and can be covered by 
the global coordinates r e R, 9 e [0, 2ir] with the periodic identification 9 ~ 9 + 2ir. 

In global coordinates the mode expansion has a discrete spectrum, 

oo oo 

n=0 l=— oo 

oo oo 

= £ £ b hlihni e- ihn ' lW - fhn ' liB + h.c. (2.10) 

71=0 1 = — OO 

where we used the null coordinates w = t + 9, w = r — 9. The weights h n j, h n j are 
related to the mode frequencies u n j and the angular momenta I by 

h n ,i = 2^ Un ' 1 ~ 

Kl = \^n,l + l) (2.11) 

and the frequencies u nt i are given by 

v nil = l + z/+ |/| +2n (2.12) 
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and v = y/1 + m 2 A 2 . 

The global vacuum is given by 

6SS|0>, = (2.13) 

and the "one particle states" are given by 

M=b%>\0) g (2.14) 

The latter transform according to the (elliptic) representation of SL(2,R) in the global 
coordinates @ 0, |23|, |£IJ. 

The Poincare coordinates cover only a finite diamond shaped patch w G [— 7T, 7r], 
G [— 7r,7r] of the infinite cylindrical boundary of CAdS 3 . The mode spectrum is 
continuous, and the operator expansion takes the form 

0= du dk b { p k e- iut+ikx + h.c. 

JO J-oo 

dh J dh b^)e- ihz - ih ~ z + h.c. (2.15) 

in the null coordinates z — t + x, z = t — x. The weights (h, h) are related to the 
energy u and the momentum k by 

h = -{uj — k) 

h = -(u + k) (2.16) 

The Poincare vacuum is given by 

O)^ ( 2 - 17 ) 
and the "one particle states" are given by 

\«>,k) = tf$\Q) p . (2.18) 
They transform according to the (parabolic) representation of SL(2,R) in the Poincare 



coordinates 29, BOL 23, 21, Bl. 



2.3 The Relation of the Vacua 

One might worry that since the global vacuum is defined over the whole infinite 
cylinder, a transformation to a single finite Poincare patch will create a mixed state. 

However, it turns out that if the parameter v is an integer, the situation is 
simpler. One can check that in fact in this case the operator O becomes periodic 
or antiperiodic over a fundamental domain which is exactly the Poincare diamond 
with periodic identifications w ~ w + 2tt, w ~ w + 2tt. In other words, we can view 
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the operator O as being defined on a torus. If the parameter v is an odd integer, 
the weights h, h are integer valued, and the operator O is periodic over the torus. If 
v is an even integer, the weights h, h are half-integer valued and the operator O is 
antiperiodic over the torus. 

The Poincare patch does not quite cover a torus (one would need more than 
one patch). To manage with a single patch, we need to cut the torus open with 
infinitesimal cuts. If we begin with a global vacuum on the torus, the cuts in principle 
will induce a logarithmically divergent entropy due to the correlations lost across the 
cut |33|, [34| . However, this contribution can be ignored. Thus, after the cuts the 



system is still in a global vacuum state. 

Now we can simply study how the modes transform under the coordinate trans- 
formation 

z = Atan(>/2) , z = Atan(w/2) (2.19) 

from the global coordinates to the Poincare coordinates 1 . One can easily check that 
the Poincare modes are analytic and bounded in the lower complex global time half- 
plane, so they correspond to positive energy modes in global coordinates 2 . 

As an aside, we also note that the transformation from global to Poincare coordi- 
nates is different from the transformation from Minkowski to Rindler coordinates in 
Minkowski space. The r = global time slice coincides with the t = Poincare time 
slice, whereas a t = Minkowski time slice is divided between two Rindler wedges 3 . 



There is one subtlety which arises from mapping the finite region in global co- 
ordinates to an infinite region in Poincare coordinates. The vacuum will develop 
a finite energy density. Subsequently the whole energy spectrum is shifted by the 
vacuum energy. We discuss this in more detail below. 

2.4 The Virasoro Generators 

The discussion in the previous section was based on treating the operators as free 
fields. We will now study the vacua again in a more general framework and examine 
how Virasoro generators (corresponding to the generators of asymptotic isometries 
of AdS3 f36|) are related in the two coordinate systems. 

The Virasoro generators are a special case of conserved charges L f corresponding 
to currents constructed from the energy momentum tensor: 

L f = J dzT(z)f(z) (2.20) 

1 Incidentally, the same transformation is used to map the infinite Minkowski space to a finite 
Penrose diagram in 1+1 dimensions (see ]35|). Thus the global coordinates are analogous to the 
Penrose coordinates and the Poincare coordinates are analogous to the Minkowski coordinates. 

2 We thank Jorma Louko for emphasizing this to us. 

3 We thank I. Bengtsson and P. Kraus for pointing this out to us. 
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where the integration is taken along an appropriate contour and f(z) is a holomorphic 
function (see e.g. paper IV.2 in j37J for more discussion 4 ). The charges satisfy the 
commutation relations 

i[L f ,L g ] = L h -~Jdzf'(z)g(z) (2.21) 

where h = f'g — g'f and the last term is an additional c-number. In the radial quan- 
tization, one takes the integration contour to be a circle of fixed radius around the 
origin, then the choice / = —iz 1 ~ n ,g = —iz l ~ m reproduces the Virasoro generators 
and the commutation relations of the Virasoro algebra with the central extension 
term where a is related to the central charge c. 

In our case, for the global coordinates on the cylindrical boundary we choose to 
define the Virasoro generators as follows: 

= f°° dw e inw T(w) , (2.22) 



similarly for the other null coordinate id we define L%>. Recall that the mode spec- 
trum is discrete in the global coordinates. 

In the Poincare coordinates, the spectrum is continuous, so we define a continuum 
of charges (similarly L^) by 

L[ P) = f°° dz e iXz T{z). (2.23) 

J — oo 

Note that L^} + 4=o * s ^ ne Hamiltonian in the Poincare coordinates 5 . Then, the 
commutation relations ( |2.21| ) imply that the L^ charges act as raising and lowering 



operators on energy eigenstates: 

[4 P) ,4 P) ] = -A4 P) - (2-24) 

Under the coordinate transformation from the global to Poincare coordinates 
w — > z(w), the stress tensor transforms as a dimension 2 operator 

T P (z) = (^\ T»{w(z)) + ^S(w, z). (2.25) 

The first term in the RHS is the classical transformation, which at the quantum 
level is supplemented by the second term which arises from the conformal anomaly. 
S(w, z) denotes the Schwarzian derivative. 



4 In what follows, we will be somewhat cavalier and not address issues of convergence in the 
definition of the charges. 

5 Compare this with the radial quantization on the Poincare plane, where the Hamiltonian is 
L_i +Z_i, see e.g. & 
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Next, we work out the Schwarzian and the transformation rules for the individual 
charges L which we defined above. The results for the latter turn out to be: 

4=0 = 4=0 + ^(4=— 1 + 4=1 ) — — (2.26) 



4Vo = E c n (A)4 9) (2-27) 

n>-l 

4< ) o = E 5 n(A)4 9) (2-28) 

n<l 

where c n (A), c n (A) are certain coefficients which can be expressed in terms of hyper- 
geometric functions of A. Note that the first relation (without the c-number term) is 
equal to the relation between the time translation generating vector field id z +idz and 
the vector field representation of the SL(2,R) generators in the global coordinates. 

Acting on the SL(2,R) invariant global vacuum |0) 5 with the 4 , we obtain the 
following results: 

L[%\0) 9 = ~ (2.29) 

4> ) olO) 9 = (2.30) 
L[%\0) g = E £n(A)^|0), ? 0. (2.31) 

n<-2 



The second equation above shows that |0) 9 is the lowest energy state, because it 

(p) 

is annihilated by all the lowering operators £\ >0 - The first equation tells that in 
the Poincare coordinates the global vacuum state is seen to have a negative energy 
density. This can be understood as a shift of the vacuum energy density associated 
with the stretching of the finite area of the patch to an infinite area in the Poincare 
coordinates. Consequently, the whole energy spectrum is shifted by the constant 
vacuum energy. This is a bit reminiscent of the Casimir effect, however in that 
case the restriction to the finite domain is obtained by imposing Dirichlet boundary 
conditions at the plates and the Casimir energy density depends on their separation. 
The last equation shows that L^ <Q acts as a creation operator on the global vacuum. 

From the bulk point of view we can formulate the result in the following way. 
It is easy to see that an observer at rest in Poincare coordinates is accelerating with 
respect to an inertial observer using the global vacuum. This leads to the naive 
expectation that the Poincare observer should see a thermal background, but as we 
have argued this is not the case. The background is not thermal but instead given by 
a shifted vacuum energy density. The acceleration of a Poincare observer represents, 
however, a critical value such that any larger acceleration will result in a temperature 



as is the case for the vacuum defined by BTZ coordinates. See []38| for related results. 

To summarize, the Poincare vacuum state can be identified with the global vac- 
uum state after the negative vacuum energy has been subtracted off. Hence in 
investigating quantum fields in AdS space, it does not make much of a difference 
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whether one starts with a global or a Poincare vacuum state. In particular, this 
result gives a more solid footing to the investigations of AdS black hole thermody- 
namics M, |2l |25|, EM . 



3. Bulk to bulk propagators in AdS space 

We now turn our attention to various propagators in the AdS space, these will be 
used later for the boundary theory description of different bulk probes. 

The propagator for a scalar field is defined as a solution to the equation 

~^d, (^ v d v G(x,x')) - m 2 G(x,x') = ^-5{x - x'), (3.1) 

where g^ v is the metric of the space and m the mass of the field. Also, x' is the 
position of the source and x the point where the field is measured. In the case of 
AdSd+i space with Euclidean signature the metric is in Poincare coordinates 

ds 2 = \ (dxl + dx\ + ... dxfj . (3.2) 

Xr • ' 



In AdS, m can be negative. A more useful parameter is given by: 



v =] J^ + mi (3.3) 

which must be real to ensure stability, this imposes a lower bound m 2 > —d 2 /4 for 
the mass |39| , 40 , f4lf . The propagator is a function of the invariant distance d(x, x') 
in AdS space, which in the Poincare coordinates is defined by 

cosh(d(x, x')) = l+ ( x o-x' ) 2 + ... + (x d -x' d ) 2 

2x x 

The solution of ( p.l|) can be found in |^2| but we include it for completeness 6 : 

. m-l 



d = 2m G u 0m (u) 



i 



T 2m\ u ') in \ 27rsinh(u) du J sinh(u) 

(3.5) 

d = 2m + 1 G» m+1 (u) = -£ Y Q v . h (cosh(tt)) 

where u = d(x,x') and Q is a Legendre function of the second kind fli3|| . These 
expressions are valid in any coordinate system provided we substitute the appropriate 
expression for the invariant distance u. 

In the Minkowski case there are different propagators depending on the time 
boundary conditions, the most relevant ones are the Feynman and the retarded 



3 This expressions are equivalent to the one in |42| which uses an hypergeometric function. 
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propagators. The Feynman propagator is the analytic continuation of the Euclidean 
one so it can be obtained from ( |3.5| ). The retarded propagator can be very different 
as is illustrated by the well-known example of 3 + 1 massless propagators in flat space. 
Recall also that in flat space the massless propagators are concentrated in the light 
cone in even dimensions. The Poincare metric in Minkowski signature is given by: 

ds 2 = — (dxl — dt 2 + dx\ + . . . dx 2 d _^\ . (3.6) 



The propagator is a function of x' , x , t—t' and r = y (xi — x'i) 2 + . . . + (x^-i — x 'd-\) 2 - 
It satisfies the equation: 

dmG+ ^ doG - daG+drrG+ d -Zl drG+ % G = l(^ 4 - Az -<W-mr) 
x r Xq 2-k~ r 

(3.7) 

In this equation the dimension d is just a parameter. A simple computation reveals 
a relation for propagators in different dimensions but for a fixed value of v: 



The only subtlety in deriving this relation is that in the right hand side of ( |3.7| ) the 
radial delta functions have different measures in different dimensions. Therefore, one 
must use the relation r5'(r) = —8(r) for delta functions acting on functions regular 
at r = 0. 

From (|3.8| ) it follows that it is only necessary to obtain the propagators in 
AdS 2 +i and AdS 3+ i. The retarded propagator vanishes for t < t' and for t > if it 
must be a linear combination of the solutions of the homogeneous equation. Among 
these solutions, we will only consider the normalizable ones. This amounts to a choice 
of a boundary condition in the boundary of AdS space which seems to be the most 
natural one given the formulation of the AdS/CFT correspondence in Lorentzian 



signature [21]. The normalizable modes can be written in terms of Bessel functions 

as 

$ = e- tujt+tkx x Jv(qx ); w 2 = q 2 + k 2 . (3.9) 



The linear combination is then chosen in such a way that 



r 



d-l 



jAr) 



- = --A^.xg-^(xo - x' )^ (3.10) 

2tt 2 r 

Using the completeness relation of Bessel functions 

I" dqqJ u (qx')J u (qx) = ^ ~ X) (3.11) 
Jo x 
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the propagator for t > t' can be written as 



G(x Q ,x' Q ,t-t',r) 



x x' f°° f°° J1 j sm(y/q 2 + k 2 (t - f)) 



IT JO JO 



dk dqq- 



cos(kr)J v (qx' )J v (qx Q ) 
(3.12) 



Performing the k integral we obtain |g3 
G(x ,x' ,t,r) = -Q(t-r)°^- 



dqqJv(qx' )J u (qx )Jo(qVt 2 - r 2 ) (3.13) 



where Q(x) is the Heaviside function (Q(x > 0) = 1 and Q(x < 0) = 0). The last 
integral can also be computed |43| in terms of associated Legendre functions of the 
first and second kind (Pg,Q%), with the result: 

1 ry\ ( ( \\ - f ( (X ~x' ) 2 + r 2 -t 2 <0, 
— P 2 _1(C0S(W)) if \ )_ , J 2 



2V2 



T y/sin(v) v ~2 



(x + x' Q Y + r 2 - r > 



G(x , x'q, t,r) = { j sin(^)_ g| i ( cosh ( M )) if / ( x o - ^o) 2 + r 2 - t 2 < 0, 



V^vrS yj sinh(w) 



(x + x' ) 2 + r 2 -t 2 <0 



if (x - Xq) 2 + r 2 - t 2 > 



where we defined 



cos(v ) 
cosh(w) 



t 2 -r 2 -{x -x' o ) 2 
2xox' 

2 -(x - 

2xqx' 



(3.14) 



(3.15) 



Note that v = d(x, x') is the timelike distance between x and x' measured in AdS 
space. 

The case when v is integer is particularly simple since sin(z/7r) = 0. In that case 
we obtain 



Gj=2 i x 0) x o> t, r) 



cos(uv) 



2tt sin(i)) 

otherwise 



if — 1 < cos(-u) < 1, i.e. v is real 



(3.16) 



To obtain the propagator in higher dimensions we apply the identity (|3.8| ). Since 
the propagator depends only on the invariant distance v (as expected) it is better to 
rewrite that equation as 

1 d 



G u d+2 (v) 



2tt sin(w) dv 



G v d {v) 



or 



d= 



d 



m— 1 



■2m\ 



(3.17) 



(3.18) 



\2ti sin(f ) dv J 

with G v 2 given by ( |3.16|) for v integer. While taking the derivatives one must be 
careful since the condition — 1 < cos(t> ) < 1 implies a discontinuity for the function 
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at cos(f) = ±1, this gives rise to delta functions. To see this better we can rewrite 
the propagator for v G Z as 



with 



Gf = L(v) 



KM = e(M 



1 

2^ 



1 d " 
2w dr\ 



m—l 



KM 



cos(n arccos(??)) 
1) = B(|r/| 



?7=cos(u) 

Tn(v) 



1 



(3.19) 



(3.20) 



where G is the Heaviside function and T n is a Chebyshev polynomial of order n. 

In odd dimensional Minkowski space the retarded propagator fills the interior 
of the forward light cone even when the field is massless. In AdS there is a further 
subtlety, as seen in fig.|I| there is a region in the future of the source where the 
propagator vanishes. Furthermore, one can show that the filled lightcone touches the 
boundary along a hyperbola that approaches the (unfilled) lightcone in the boundary 
of even dimensional Minkowski space as the source approaches the boundary. 

The calculation proceeds in a similar way in the case of 
AdS 3+ i MM. The expansion in modes reads: 




G(x ,x' ,t,r) 



[x x ) 



7T 



JO 



dkdqkq ^-^=^ -^- J ( Ax) J„ ( qx' )J u (qx ) 
q (3.21) 



Performing the k integral we obtain M3 

G(x ,x' ,t,r) 



-G(t-r) 



{ x o x o) 
2tt 



o 



dqqJ u (qx' )J u (qx ) 



cos 



(3.22) 



Figure 1: For d = 2n + 
1 and v £ Z the prop- 
agator is non-vanishing 
only in the shaded re- 
gion. Here 77 = cos v, see 
eq-CT). 



This integral is divergent. To handle it we can rewrite it as 
^IjLj^dq sm(qOJv(qx' )J u (qx ) 



X0X'„ 



2tt (d( 



ro 

cos(f7r) 



Q, 



x o+ x 'o £ 2 
2xox' 



£=Vt 2 -r 2 

if £ < \x' - x \ 

if \xq — Xq J <C £ <C Xq ~\~ Xq 

if Xq + Xq < £ 

(3.23) 



In this case the propagator is simpler for half-integer v 
(which includes the case m = 0). Using the notation of 
eq.( |3.15p the propagator can be written 



G veZ+l Hxn x' t r) 1 dKM 



r]=cos(v) 



with 



KM = e(M 



l )~ P v 1 



(v) 



(3.24) 



(3.25) 
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Again, using ( |3.8| ) we can express the propagator in higher dimensions as: 

(3.26) 

r]=cos(v) 

Moreover, in the case we are looking at, namely v — ~ integer, the function P v _i is 
a (Legendre) polynomial of order v — | and so the derivatives terminate for a given 
d. This means that the propagator is concentrated on the light-cone, since the only 
derivatives that remain are due to the jump of the function. 

Finally, all the propagators are independent of a choice of a coordinate system. 
Since they transform as scalars, one only needs to express the invariant distance in 
the desired coordinates. 



4. Examples 
4.1 Bubbleography 

Let us now use the propagator which was derived in the last section to study the 
boundary theory description of prototypical probes of the bulk geometry. As a 
concrete example, we consider a massive source falling along a geodesic in the AdS^+i 
space. In global coordinates the metric is given by 

ds 2 = — cosh 2 /jdt 2 + d/j 2 + sinh 2 n,dVt 2 d _ l (4-1) 

and we will consider a source moving along a geodesic at // = 0. Primed coordinates 
will refer to the source, while unprimed coordinates will refer to the point where we 
evaluate the field. We are interested in the leading behavior of a field that couples to 
the source close to the boundary where — > oo . To find this we need to integrate 
the propagator along the worldline of the source, evaluate it in a point a little bit off 
the boundary and take the limit /i — > oo and extract the relevant term. The integral 
will get contributions from a small segment of the worldline, < v < 7T and hence 
we find 

<f> ) = -— / dsl—^—) , (4.2) 

Ltx J \2tt sin v dv I sin v 

where n = d/2. We will only consider integer values of v. For the // = worldline 
ds = dt', and the invariant length is given by 

cos v = cosh fi cos (t — t') . (4.3) 

This implies 

dt' sint> 1 , 

— — — (4.4) 

dv cosh \x R — cos 2 

V cosh /i 
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and hence 



, . 1 f n 7 smv I I d\ cosuv , A . 

4>(Q,fj) = —— dv - — : — . (4.5) 

2ir Jo cosh u Ji— cos t v V 2tt sm v dv I smv 

V cosh 2 /i v ' 

The integral over v will pick out powers of cosh \x such that the leading contribution 
for /i — > oo is 

= A ni! ,(cosh/i) _n_1/ , (4.6) 

where A njV is a constant that is zero for v + n — 1 odd. Hence we conclude that the 
boundary field is constant. 

In the Poincare coordinates the situation is different. The worldline is now given 

by 

x' 2 = l + t' 2 , (4.7) 

implying the line element 
We find 

1 r n , I dt' ( 1 d \ n ~ cos z/-u 
4>{x,x Q ) = -— dv- -— - — : — . 4.9 

If we use the expression for the worldline in the expression for the invariant length 

cos^^fa-^-C-'') 2 (410) 
(we have chosen x' — 0) and define 

*' 

sin a = . (4.11) 

sin/3 = , * X ° X (4.12) 
+ (t 2 - x 2 - x 2 - l) 2 

we obtain 

— 2x cosv 

sin (« + /?)= (4.13) 
+ {t 2 -x\-x 2 - if 

We can now obtain (noting that (5 is independent of v) 

1 r* da ( 1 dY^cosvv 

cf)(x,x ) = - 7r / ^— — — — (4.14) 

27r Jo ctf \2,ii smv dv J smv 

f n , a ■ 1 ( 1 d \ a l cosuv , A . 

= -/ dvAsmv—7 I (4.15) 

Jo ^1 — A 2 cos 2 v \2n smv dv J smv 
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where 



A= , " 2 "° (4.16) 



At 2 + (t 2 -xl -x 2 - 1) 



The v integration results in 



B n>u ^ _ (4.17) 



which implies a boundary field expectation value of the form 

(O) = C n , v ^, (4.18) 

(At 2 + (t 2 - xl - x 2 - l) 2 ) 2 

where 5 n ^ and C n ^are constants that are zero for v + n — 1 odd. As an example of 



the scale/distance duality between the boundary and the bulk, in [|22| it was shown 
that an object at a fixed distance from a horizon would look like an extended blob 
in the boundary theory; the closer to the horizon the larger it will appear. A freely 
falling object is then expected to grow as it closes in on the horizon. The above 
calculation shows that this is not the whole story. Even if the boundary description 
of the object starts out like a blob it will develop into an expanding bubble over 
time as illustrated in fig.|2]. This can be understood intuitively in the following way. 
Since the object is falling, all scales will expand but at different rates. The points in 
the center of the bubble depend on the source at later times, when the information 
has not had time to spread much in the transverse directions. On the other hand, 
the source is falling faster at later times and therefore corresponds to a region in 
the boundary that must expand faster. The result is a shockwave that produces a 
bubble. 

In the case of AdS3 we can also verify that the change of coordinates from global 
to Poincare amounts to a conformal transformation in the boundary; the bubble is 
hence a conformal transformation of the constant profile. In fact, 

At 2 + (t 2 - x 2 - if = (1 + z 2 )(l + z 2 ) (4.19) 

(with z as in the previous section) and since 

dw 2 



dz 1 + z 2 

we find that the bubble is simply given by the conformal scaling prefactors 



(4.20) 



n-\-v 

(1 dw dw\ 2 
4&a?) ' < 4 ' 21 > 

As we see, in the case of AdS^, a coordinate transformation in the bulk induces 
a conformal transformation in the boundary and the boundary field transforms ac- 
cordingly. The geodesic // = where the particle is always at the center in global 
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coordinates was particularly simple. However any other geodesic can be obtained 
from that one using an isometry of AdS^ space and so the boundary field produced 
by such particle can be obtained by the corresponding conformal transformation in 
the boundary. 

In Poincare coordinates the trajectory is x'q = 1 + t' 2 . Performing the isometry 
— *■ x /a, t — > t/a,x — > x/a, with a an arbitrary parameter we obtain the new 
trajectory defined by 

x' 2 = a 2 + t' 2 (4.22) 

In the boundary, the conformal transformation induced by these coordinates is simply 
a rescaling of z, z. However in global coordinates the conformal transformation is 
not so trivial and is given by 



w'\ / w\ (w'\ I w 



tan y~2 J = fl tan y~7^j > tan ^— j = a tan y^j (4.23) 

where vJ are the new coordinates. Since in the coordinates w the boundary field is a 
constant we obtain the field produced by the particle following the new geodesic as 

(4-24) 



dw' I \ dw' 



a? +1 



(a 2 cos 2 f + sin 2 f ) 2 (a 2 cos 2 f + sin 2 f 



(4.25) 



As a check it can be seen that the same expression is obtained if one plugs in the 
new geodesic in the previous formulas with the propagator. The resulting field looks 
for small times also like an expanding bubble but is of very different nature since it 
is obviously periodic in time. 

A more interesting case is that of BTZ coordinates. The conformal transforma- 
tion from the coordinates w where the field is constant to BTZ coordinates reads: 

tanh ^(r + — r -)^7^j = fl tan > tanh ^(r + + r -)~7^ = °tan ^— ^ (4.26) 

and the field is accordingly 

a u+u r 2 _ r 2_y-±± 



2 



(a 2 + (1 + a 2 ) sinh 2 2 ( G 2 + (1 + a 2 ) sinh 2 

(4.27) 

In the black hole case, we must perform the identification <ft = <ft + 2n which amounts 
to a sum over images: 

oo 

BTZ (0 + ,0_)= <K0+ + 27m,0_+27m) (4.28) 

n=— oo 

The field is depicted in fig.|| where <fi = —it . . . tt. 
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Figure 2: Dilaton field in the boundary of, left: AdS% for v = 2 and times t = 1,2,3 and 
right: a BTZ black hole for v = 2, r+ = 1, r_ = and times t = 0.5, 0.8, 1.2. 



4.2 Supersymmetric YM at zero temperature 

From equation ( f4.6|) it follows that the above picture is no longer relevant if v — n + 1 
is odd. In these cases the expectation value in the boundary theory vanishes. In 
particular this is true for massless dilatons, which in the case of AdS§ implies that 
the expectation value < F 2 > must vanish for a bulk source in a free fall. How can 
one then understand the nonvanishing expectation value obtained in [2!2]j ? 

The important point is that in \2~2\ the source is held at a fixed position in 
Poincare coordinates. This can be achieved by suspending it on a string hanging 
down from the boundary with an appropriate tension. The other endpoint of the 
string will then show up as a point charge in the boundary theory. We conclude that 
nonvanishing expectation values of F 2 must be supported by charges. To be more 
precise we can use the propagator derived above to calculate (F 2 ) due to a point 
charge by representing the charge by a string that goes straight down to the horizon. 
We must then consider the string as an extended source for the dilaton field and 
integrate the propagator along the string. 

To perform our calculations we must know how F 2 couples to the supergravity 
fields. There are two such fields that will be of interest to us, the dilaton <fr, and the 
volume. We will be interested in variations with respect to certain linear combina- 
tions of these fields and a convenient basis is to consider variations either with the 
string metric or the Einstein metric held fixed. To fix the conventions, let us consider 
the D-brane action 

I BI = TrT p J d p+1 xe-^det(g^ + T^ v + ij'rO.X ,<) V X ;,). (4.29) 

Following e.g. 0] we then rescale the worldvolume gauge field strength, 

T = TF , (4.30) 

where T is the string tension, 
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and substitute the Dp-brane tension 

T p = g; L (2n)^T^. (4.32) 
Expanding to second order in F we find 

I BI = Tr [ dP +1 x^e-* (t p + -\- {g^F^f + ... 

J \ 4g Y M , 



(4.33) 



where the Yang-Mills coupling is given by 

9ym = {2*y- 2 g s a' v -^ . (4.34) 
It is convenient to define a new metric g^ u , defined through 

9»v = e^g^, (4.35) 



which for p = 3 coincides with the Einstein metric. The Born-Infeld action then 
becomes 

I BI = Tr f d p+1 xJ^d (t p + -r-^-e"^ {g^FJf + ■■) ■ (4-36) 



4# 



YM 



If we furthermore define 

J=— (4.37) 



we find that F 2 is obtained as 

5hl ~ 1 -F 2 . (j.38) 



50 4^ M 

For p = 3 and p = 4, we have checked that <fi fluctuations keeping g^ v fixed are 
massless. On the supergravity side, the supergravity action (dimensionally reduced 
over the unit sphere) is then given by 



S = -T 5 ? / cF +2 xy/^g%dJd u l (4.39) 

4/C 10 J 



where Qs-p is the volume of the unit 8 — p sphere. If we have a string as a source we 
will need to vary 

s - = db / = / ^ ^ (4 ' 40) 

where the metrics are now the induced ones on the worldsheet. We see that the 
variation with respect to gives 



5S„ 



hit 
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The corresponding variation of the supergravity action is then given by the boundary 
term 



5S 



8 ~p / jp+i 



2k 2 



10 



dP^x^gg uv eT^d { 



v 



(4.42) 



where is the background value of the dilaton. Hence it follows that 



4<7 



YM 



,2\ _ 

2k w 



(4.43) 



17->oo 



if we consider string like sources. 

Let us apply this to the case of a string hanging straight down when p = 3, i.e, 
d = 4. The string traces out a worldsheet that gives rise to a field 



<f>(U,x) 



1 2k 



10 



47ra' tl 5 

If we want to write the AdSs metric as 



dt'dU'yf=g~ 2 



d cos 2v 



2ir sin v dv sin t> 



(4.44) 



ds 2 = a' ( (J) 2 df/ 2 + (£) 2 ^ dz 2 - dt 2 ) + R 2 dn 



(4.45) 



where i? 2 = Qym^/N we need to take x — > x/R 2 , t — > t/i? 2 , rescale the metric by 
a'^ 2 and put x = 1/E/. Noting that 5 (3) (ar) = ^(J® (x/i? 2 ) we find 



1 2k 



10 



dvdU' 



d cos 2t> 



where we have put r 2 = I^f =1 x 2 and also used 



[7-2 + jji-2 + 2U~ 1 U>- 1 cosv 



t-t' 
R 2 



ijj-2 + v ,-2 + — - 2U-W- 1 cosv 



dv smv 

(4.46) 



(4.47) 



We find 

4>(U,x) 

hence 



15 2^, u-<u-< 



256vr 2 J R 6 fi, 



(u>- 2 + £ 



r 2 \7/2 



' :R 2 U- 41 



1287T 2 



4 ' 



4t? 



KM 



g YM VN^ 



32tt 



(4.48) 



(4.49) 



which suggests that the effective charge goes like ygyMvN . Hence the force on 



another charge y gyM^/N is ~ gyM\fN \ implying a potential ~ gyM^N-, in agree- 
ment with [0. 
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4.3 Supersymmetric YM at nonzero temperature 

We may also consider the nonzero temperature case. The starting point is the metric 
of a non-extremal p-brane, given by 



ds 2 {R\*dU 2 o {U\^ " 

f(U) = 1 - ^ (4-50) 



while the dilaton is given by 



We have defined 



3-P 



= 9s*'- ( 4 (4.51) 



R— = 9ym \j dpN (4.52) 

where 



dp = 2 7 - 2 %^r(^) (4.53) 



9-3 P t _ , 7 — p, 

2 

and gyM is given as before. 

We will limit ourselves to the case of d = 4 and begin by deriving the propagator 
for a static configuration since we can no longer use the exact retarded propagator 
evaluated in the previous section. 

The metric ( 4.50|) reduces to 

ds 2 = jj^dU* - ^j^dt 2 + ^dx 2 + R 2 dQl (4.54) 

f{U) = l-{^) (4.55) 

In this case we will consider only the case of a static source coupled to a massless 
field through the action 

5 int = / ds(j){x{t)) = I dty/\^\<f>(x{t)) (4.56) 



In the propagator equation ( |3.1| ), this amounts to multiplying the right hand side by 
g tt \ and omitting the time derivatives. The resulting equation is 



(4.57) 

The boundary conditions are that G is regular at U = Ut and that G vanishes at 
U — > oo. Rescaling the coordinates as U — > UtU, t — > tR 2 /Ut and x — > xR 2 /Ut, the 
equation simplifies to 



(u 5 - u)d uu G + (5m 4 - l)d u G + ud u G = 1 S(u - u')S {3) (xA (4.58) 

R'u 
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Fourier transforming both sides of eq. ([4.57| ) with respect to Xi leads to 



~ \/?y 4 — 1 

(u 5 - u)d uu G + (5m 4 - l)d u G - uk 2 G = ^— S(u - v!) (4.59) 

R'u 

with 

G(u, u', k)= [ d 3 ke l ^G{u, u', z<) (4.60) 



For u < u' and for u > v! the Green function is given by (different) solutions of the 
homogeneous equation 

(w 5 - u)y"{u) + (5m 4 - l)y\u) - uk 2 y(u) = (4.61) 



The homogenous equation has been obtained in [|17j] and discussed at length in [|46|]. 
For u > u' there are two solutions which we shall call yi(u,k) and y\(u,k) which 
satisfy 

y 1 (u,k) = \ + 0(l E ), u^oo (4.62) 

= l + 0(±^), u^oo (4.63) 
u u 

(4.64) 

and for u < u' there is another pair of solutions 2/2 (w, k), jj2{u, k) satisfying 

y 2 (u,k) = l + 0(u-l), u^l (4.65) 
y 2 {u,k) = ln{u-l) + 0{u-l,{u-l)ln{u-l)), u^l (4.66) 

(4.67) 

The boundary conditions then require that 

G(u> u',k) = A yi (u,k), G(u<u',k) = By 2 (u,k), (4.68) 

and A and B should be chosen such that 

By 2 {u',k)-Ay 1 {u , ,k) =0 (4.69) 
By' 2 (u', k) - Ay[(u\ k) = )== (4.70) 



It follows that 



^4 = y i 1 C4 71) 

yiy'2-y2y'iR 7 u 2 V^i 

B = T ] (4.72) 

2/12/2 - 2/22/i i?^ 2 ^^? 3 ! 
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where the functions are evaluated at u = v! . As follows from ( |4.61| ), the Wronskian 

w(y\, y 2 ) = 2/12/2 - y^y'i is s iven b y 

W(y u y 2 ) = ^- (4.73) 
u — it 

where w(k) is a function that depends only on k. This simplifies the expression for 
the propagator which turns out to be 



1 vV 4 - 1 

G(u, u, k) = ^j^yi( u >i k )V2{u<, k) ^ /R7 (4.74) 

with m< = m.m(u,u') and w> = max(«,M'). In ref. [^E] a series expansion was 
obtained for yi and y 2 and also for w(k). There it was shown that for real k there 
are no zeros of w (and so the propagator is well-defined) and also that there is an 
infinite number of zeros on the imaginary axes. These zeros give precisely the masses 



of the glueballs in [17 



1 f 1 f UjR 

S int = ^ J dudt V9ttguu<j> = J dudt^ t —==(p (4.75) 



For what follows, besides the expansion of [46| it will be useful to have also an 
approximate but simpler expression for w(k) which we derive in the Appendix using 
the WKB method. 

A charge in the boundary is represented in the bulk by a string hanging from 
the boundary down to the horizon. The minimal action configuration corresponds 
to the string extended only in the U and t directions. Such a string will be a source 
for the dilaton through the coupling 

1 f .„ ,. 1 /"...__ uR 

where we have put a' = 1 for simplicity. The boundary value of corresponds to 
the expectation value of F 2 in the presence of the charge. We know that for short 
distances it will behave as 1/r 4 since in that region we can replace the black hole 
backround by AdS space and so ( |4.49| ) is reproduced. For large distances it should 
decay exponentially as exp(— r/Ls) where Ls is the screening length. The value of 
the dilaton can be computed using the static Green function as 

*<) = ±J du J ^e-^-0=G{u, u', k) (4.76) 
Using ( [4.74 ) the value near the boundary is given by 

<p(u,x)^oo ^-M X ) ( 4 - 77 ) 

' AnR^J J (2vr)3 e w(k) 1 ' 

The angular integral in the Fourier transform can be performed with the result 

Mx) = j2^J dU l dk kr w(k) (4 ' 79) 
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where r 2 + x\ + x\ + x\. The integral in k can be extended from — oo — > oo and 
evaluated by the method of residues. The integrand has simple poles in the zeros 
of w(k) which are located on the imaginary axis. A zero in w(k) means that the 
solutions yi and y<i are proportional, that is that for that value of k there exists a 
solution well behaved in u — 1 and u —>■ oo. Such solutions are the eigenfunctions 
computed in [|17|, [46| and so, the zeros of w{m n ) = define the glueball masses in 



Jj4 oo g-m n r 

MX) = 8^ £ mn rw>(m n ) J dU ' V2{u '' mn)dU ' (4 ' 80) 



the 2 + 1 confining theory. 

E#_ 

where m n and y2{u',m n ) are the rz eigenvalue and eigenfunction of eq.( |4.61| ). Note 
that the integral in u' is well defined for the eigenfunctions. Approximate values of 
the eigenvalues and eigenfunctions are obtained in the Appendix. 

At long distances the function decays exponentially with a screening length given 

by 

Ls = j^— (4.81) 
U T m gL 

where m g \ ^ 3.4 is the lowest eigenvalue. From the boundary point of view this was 
expected, since at large temperatures the theory is reduced to one less dimension. 
The mass gap in this dimensionally reduced theory is m g i which in the original theory 
is interpreted as a screening mass. 

Hence we find that the screening length is determined by the glueball masses 
and is of order 1/T = R 2 /Ut- 

4.4 Confinement and flux tubes in non supersymmetric YM 

The metric fl4.50| ) can also be used to describe nonsupersymmetric QCD P , fI5| . To 



achieve this in the metric ( |4.50|) we must go first to Euclidean time t — > r and then go 
back to Minkowski signature but through x p — > t. To avoid a conical singularity at 
U = Ut, t must be periodic with period tt/T = ttR 2 /Ut ■ As a result the boundary 
theory is a SYM theory compactified on a circle and with coupling constant gy M T, 
where gy M is the coupling constant of the higher dimensional theory. 

The potential between two quarks in the boundary is described by the minimal 
action of a string going from one quark to another and extending in the bulk. When 
the distance between the two charges is very large, the string hangs down to U = 
Ut where the metric is approximately flat and the energy is proportional to the 
separation length giving confinement. This string is a source for the dilaton and so 
in the boundary is seen as a non vanishing (F 2 ) which will be concentrated along a 
flux tube connecting the two quarks. 

The tension of the flux tube was calculated in |§ using the picture of a hanging 
string with the result 
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This can be obtained also in a different way. The free energy A(/3, A) and the gluon 
condensate are related by 

where we defined A = g\ M - In this case the gluon condensate is known from the 
boundary value of the dilaton and so the free energy can be obtained up to a constant 

as 

A = -/t/^4T (F2) - (4 ' 84) 
To proceed, we first need the dilaton equation 

' d^g^e-^^-^V^S. (4.85) 



24,"^ JJ " 4vr« 

Let us integrate over the 8 — p dimensions of the internal sphere and the p — 2 
transverse directions to the flux tube, i.e. all directions transverse to the string 
except the radial U and the angular r. We put x = I d 6 x(j) and obtain 

7-p 

1 -duV^9 UU e-^d uX = ^(^) 2 5 {2) (U - U T ) (4.86) 



2/«f 47T \R 



Integrating over the (U, r) plane then gives rise to a boundary term at U — > oo. Note 
that the origin of the radial U coordinate at U — Ut (the horizon of the Euclidean 
black hole) is a regular point since the period T has been chosen in an appropriate 
way. We find 

^ / <T-*xdTy/=gf u e-*d a $ =^{^Y ■ (4-87) 



2 K io 



Note that the leading contribution to (f) is constant over the S 8_p near the boundary. 
(In the case p=3 one has that the Kaluza-Klein states on the sphere correspond to 
states with higher conformal dimension). Using the previous results we finally find 
that the integrated gluon condensate is given by 



7-p 

1 fU T \ — 



d p - 2 xdr^^F 2 = — 9 / (P-'xF 1 = — — . (4.88) 

We must now relate the gluon condensate to the tension of the flux tube. In the 
perturbative case the gluon condensate is proportional to A and the energy and gluon 
condensate are equal. The gluon condensate calculated above is, however, given by 
-^(F 2 ) ~ y/\. The reasoning is identically the same as in the case of ( [4.401 ). The 
result is a free energy that is twice as large as the gluon condensate 7 . This implies 



7 For a D-string the factor is —2. 
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that ( |4.88| ) correctly reproduces the tension ( f4.82| ). More generally we see that, for 
massless fields, the total energy in the boundary theory is obtained by integrating 
the source over the bulk. 

In fact (F 2 ) can also be computed as a function of the position and that gives 
the profile of the string. We first will look at the case of a 2 + 1 dimensional theory 
and then at the 3 + 1 case. 

When p = 3, the metric is given by 

ds " = T7T7rFK du2 + ^P^ 2 + S(^r + d A ~ dt 2 ) + R 2 dQ 2 5 (4.89) 
U A j{U) R A R A 



m) = l -{ir) ( 4 - 9 °) 

where we have let x% become our new time coordinate t. A difference with respect to 
the non-confining case is that at U = Ut there are geodesies of constant U. In fact in 
that region the metric is approximately a product metric of an R 2+1 parameterized 
by %i,2,t, a five sphere and a plane parameterized by U,t with \JU — Ut a radial 
coordinate and r an angle. From now on we will consider fields that are independent 
of the coordinates of the sphere. 

A particle sitting at the center U = Ut of the (Euclidean) black hole will remain 
at rest. Furthermore, a freely falling source will eventually be trapped around the 
R 2+1 at the center and henceforth represent a stable configuration. This is the reason 
for the existence of a mass gap and stable glue balls as well as of confinement. If the 
particle is a source for the dilaton, then it will induce in the boundary an expectation 
value of (F 2 ) which will look like a stable blob and be interpreted by the boundary 
observer as some kind of a glueball state. Consider now the more physical case of a 
string sitting at the center of the (euclidean) black hole and extended for example in 
the x\ direction. In the boundary it will give rise to a flux tube with a shape that 
we proceed to calculate. 

First the dilaton field produced by the string should be obtained. The field 
is independent of t and X\ since the source is static and extended in x\ and is also 
independent of r since the string is at the center of the p = y/U — Ut, t plane. The 
equation Q4.85|) for (f) is: 

da (u 5 f(U)dum x 2 )) + R 4 Ud 22 m x 2 ) = ^ ^ S(U - U T )5(x 2 ) (4.91) 

where S7 5 is the volume of the unit five-sphere. Rescaling as before U — > Utu, 
t — > tR 2 /Ut and x xR 2 /Ut, the equation simplifies to 

(u 5 - + (5m 4 - l)d u( j) - uk 2 ^ = L_^ fo - 1) (4.92) 

where the Fourier transform 

cf)(u,k)= f dke ik x 2 (f)(u,x 2 ) (4.93) 
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was introduced. The homogeneous equation is the same as ( |4.61| ). In terms of its 
solution the field can be written as 



(P(u,k) = Ay 1 (u,k) (4.94) 

The solution yi(u,k) behaves as ln(ii — 1) when u — > 1 as is needed to match the 
delta function on the right hand side. In fact we need 

^)^i|f^ln( U -l) (4.95) 
On the other hand we have 

yi (u, k) = ay 2 (u, k) + /3y 2 (u, k) (4.96) 
p = y\V'2-ViV'i = W(y 1 ,y 2 ) = w(k) 
V2V2 - 2/2^2 W(y 2 , m) 4 



where w(k) is the function defined in ( [4.73 ) 



The constant A can now be obtained since 

(f)(u, k) = Ayt{u, k) = ay 2 (u, k) + fiy 2 (u, k) «w Af3 ln(w - 1) (4.98) 

= ~ 4tt U T R 2 Q 5 w(k) (4 ' 99) 
The field created by the static source is then given by 

*<"■*> = -^imrj wf ,tx W) vMk) (4100) 

Using that near the boundary y\ pa u~ A = U^/U i and using the relation ( |4.43j ) 
between the boundary field and the gluon condensate, it follows that 

1 (^H^/Ae-^^TT (4.101) 



4<# M N ' AtiR^J (2tt) w(k) 

This will be seen by the boundary observer as an infinite string of finite width. In the 
case of two charges separated by a large distance the string will join both charges. 
As a check, the integral of this profile over the transverse direction x 2 is given by the 
value of the Fourier transform at k = 0. Since w(0) =4 (as can be deduced solving 



the differential equation exactly for k — 0) the result ( 4.88 ) is reproduced. 

The profile of the QCD string is easily seen to decay exponentially at large 
distance as exp(— x 2 ) with the lowest glueball mass. An approximation to 
the profile can be obtained using the function w(k) computed in the Appendix and 



is 



4[/| r dk 71 
^9ym^ ' ~ 4:nR 4 J (2tt) l6y/2 cosh(ak) 

7T 2 m 



( F ' 2 ) « i^L I — V 77 7k J ZZTTZIR ^ (4-102) 

(4.103) 



2f RW cos h 3 (gf ) 
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where in the last line we return to the original variable x 2 —* Utx 2 /R 2 and a is the 
constant defined in the Appendix as 

1 / 1\ 2 

« = 7^r: ( 4 - 104 ) 



4V27T V 4 

This profile is depicted in fig.©. 

For the p = 4 case the metric and expectation value of the dilaton are given by 

ds2 = {vfw) + (1) 1 m)dr2 + (i) 1 ( ^ dA ~ dt2) + RiuUnl 

e* = 9s(^y (4.105) 

Considering a string situated at C/ = Ut and extending along x$,t as a source, the 
equation for the massless scalar <fi can then be written as: 

du (u 4 (l - ^pj du$\+U^du$ = ^2k 2 w V^^5 {2 \x 1:2 )5(U-U t ) (4.106) 

where is the volume of the 4-sphere. Rescaling the variables as U — Utu and 
Xi — > XiRz I ^/Ut the equation reduces to 

(u 4 - u)d uu $ + {Au 3 - l)d u 4> - uk 2 4> = ^$LJ=5{u - 1) (4.107) 



47rft 4 R 3 VTh 

where the Fourier transform 



4>(u,k) = J d 2 ke il < g $(x 1 ,x 2 ,u) (4.108) 



was introduced. The equation is similar to (|4.92|) . Two solutions of the homogeneous 
eqution can again be defined by 



yi(u, k) = u 3 + 0(u " ; // — :x 



(4.109) 



y 2 (u,k) = H-O(u-l) it — > 1 
The Wronskian between these solutions is given by 

W( yi (u),y 2 (u)) = 2/12/2 - Wi = (4.H0) 

u — u 

where w(k) is a function only of k and can be computed as a series expansion as was 
done in or in a WKB approximation as we do in the Appendix. Now the field 
can be written as 
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Expanding this field near the boundary u — > oo and using the relation ( |4.43| ) with the 
expectation value of F 2 , the Fourier transform of the string profile can be obtained 

as 



4^< F2 >^ = i(i; w {k) 



(4.112) 



Again one can check, using that w(0) = 3, the integral of this profile reproduces 
( |4.88|) . The string profile is given by 



d 2 k- 



4tt V R 



w(k) 



a 5 Jo 



sinh ny 



dy 



(4.113) 
(4.114) 



where r 2 = x\ + x\. The last integral cannot be obtained in closed form, but two 
expansions valid for large and small r are: 



(y 4 + yVo( 2 

sinh ny 



dy 



TXT 



9 oo , 

-J2(-) n n 2 (n 2 -l)K (n- 



n=2 
^ oo 



.27) 



E 



(4.115) 
(4.116) 



c n = 2(2 5+2n - l)r(5 + 2n)C(5 + 2n) + 87r 2 (2 3+2n - l)r(3 + 2n)C(3 + 2n) (4.117) 

From the first series it follows that at large distances the profile decays exponentially 
as exp(— 2-jrr/a) which is right since 2ir/a is the mass of the lowest lying glueball in 
the approximation we are considering. The profile obtained is depicted in figure ^. 
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Figure 3: Profile of the QCD2+1 (left) and QCD^ + \ (right) flux tube. The r-axis has 
units of the corresponding l/m g i and the profile is normalized so that its integral is 1. 
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5. Conclusions 



We have studied various aspects of the bulk/boundary correspondence. First, we 
addressed a basic issue about the possible vacuum choices in AdS space. Often 
one imposes the system to be in a Poincare vacuum state, this vacuum is however 
defined with respect to a coordinate system which does not cover the whole manifold. 
Hence one might wonder whether it would be more natural to use the global vacuum 
state, which is defined with respect to the global coordinates covering the whole 
space. However, we showed that for practical purposes the two vacua can in fact 
be identified after subtracting off a constant zero point energy contribution to the 
vacuum. 

We then turned our attention to various probes of the bulk and their description 
in the boundary theory. We started by an extension of an example considered in 
Iffil , a source in AdS space. Instead of the quasistatic treatment of |22| , we took the 
source to be in free fall, and evaluated the projection of the source to the boundary 
by using retarded propagators instead of static propagators. This results in a slight 
modification of the results in [22]; instead of looking like a blob with a single scale, 
its size, the boundary profile looks like an expanding bubble characterized by two 
scales which encode the radial position and the state of motion of the source falling 
in the bulk. 

Static sources are also possible, but in that case they should be kept fixed by 
e.g. suspending them on a string. The presence of the string is in turn represented 
by an additional point charge in the boundary, which induces a nonzero expectation 
value for the field strength, (F 2 ). We evaluate this expectation value and find the 
correct form for the potential from the point charge. 

At non-zero temperatures, the point charge is thermally screened. We examined 
this well, and evaluated the thermal screening length. It turned out to be 

given by glueball masses, of the order of the temperature T. 

We also studied non-supersymmetric, confining QCD. In that case, the center 
of the (Euclidean) black hole plays a central role. The real world at the boundary 
has the character of a shadow of the world at the center of the Euclidean black 
hole. A string sitting in the center, appears in the boundary as a flux tube of finite 
width with a profile that can be computed. Pushing this picture a little further, 
a closed string near the center of the black hole is mapped into a glueball in the 
boundary. If the string state is light, then the center of mass will be spread around 
U = Ut which is the case studied in Ul7j . One can also consider a qualitative picture 
of scattering between such glueballs 8 . If the bound states scatter at high energy, 
two distinct regimes appear. At small scattering angle the momentum transfer is 
small and so the particles in intermediate channels are the lowest lying glueballs. 
In that domain, when states far from the U = Ut do not contribute, one expects 



3 We thank Bo Sundborg for a discussion on this point. 
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a Regge behavior from the string scattering amplitudes which corresponds to the 
expected Regge behavior of glueball scattering in the boundary. At finite angles 
the momentum transfer is large and the region in which the bulk metric is AdS$ is 
expected to dominate. In the boundary that should correspond to the parton region 
in QCD, but unfortunately the theory studied at high energies is not QCD but 
N = 4 SYM, which is precisely what string theory in Ad S§ describes according to 
Maldacena's conjecture. 
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Appendix 

In this Appendix we will compute, using a WKB approximation 9 the function w(k) 
defined in eq. (|4.73|) . To obtain a WKB solution of eq.( |4.61| ) it is better to work 
with imaginary k and then analytically continue to the real axis. In that case, the 
equation (|4.61|) can be recast in the form of a Schrodinger equation for a zero-energy 
eigenstate fT7|| : 

- x "(x) + V(x)x = (5.1) 

. . k 2 3x^ Qx 2 1 , . 

V(x) = —-r—rr, tt+ -. 9 , 9 tt^ 5.2 

Ax(x 2 - 1) Ax 2 {x 2 - l) 2 



where x = u 2 and x( x ) = \Jx{x 2 — l)y(-\fx). There are two regions where the 
solutions of interest can be expressed in terms of Bessel functions : 

Xi(x, k) = ^x(x 2 - k) « ^ J a (^j) ifx»l 

X2 (x,k) = ^x(x 2 -l)y 2 (^,k) « p(x-l)Uk^±) if(rr-l) < 1 

In the region (x — 1) ^> 1/k 2 , x <^ k 2 /3 the potential is approximately given by 

For large k we can use a WKB approximation: 

X {x) w -^j ^e < X.°V^W + ^rv^wj (5.5) 



9 while this paper was being written ref. |47j| appear which has considerable overlap with what 
follows. 
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The strategy now is clear. For large k there is an overlap between the region x—l^> 
1/k 2 and x — 1 -C 1, so we can compute the coefficients A 2 B 2 corresponding to the 
function X2- There is also an overlap between the regions x»l and x <C k 2 /3, so 
we can compute the coefficients A±, B\ corresponding to xi- Finally, the Wronskian 
between the two WKB solutions can be computed. The coefficients are given by 



An = B* = * e ^~ ife V^+*JT\/^) 

z 2 V2tt 



(5.6) 



where we assume that the two points a, b are such that 1 ^> a — 1 ^> 1/k 2 and 
1 <C b <C k 2 /3. The Wronskian can now be obtained as 

W( X i, X2) = ^ = 2i(A x B 2 - A 2 B 1 ) (5.7) 

8^ ( k r b dx ,. , . ... 

■cos - / . +-= + k\ ——\ . (5. 



7T 



k 2 




The argument of the cosine is approximately independent of a and b, so we can take 
b = 1, a = oo. The result is 

w(k) = — ^^-cos(afc) (5.9) 



7T 



k 2 



a - 1 [°° dX - 1 (v( l ) ) r 10 1 



From here the masses of the ++ glueballs follow immediately as 

71 1 

m = -(n + -) (5.11) 
a 2 

The obtained values are in good agreement 10 with the numerical ones except for 
the fact that there is an extra eigenvalue at n — 0. Of course this is at small k, 
i.e. outside the validity of the approximation. However a simple modification of 
the Wronskian cures this problem and gives a good approximation in the full range 
k = — > oo for the function w(k) obtained numerically. The modification is to write 

w(k) — -, — — cos(afc) (5-12) 

» (* 2 - (£) ) 

In this approximation, w(0) ~ 5.017 to be compared with the exact value w(0) = 4. 
Analytically continuing k — > ik we obtain the desired function 

w(k) = — j - 2X cosh(o:fc) (5.13) 



7T 



10 As it turns out, the values of pTfl are even better since there a next-to-leading order was 
considered 
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The other case studied in the text corresponds to the homogeneous equation 

(u 4 - u)y"{u) + (4m 3 - l)y\u) - uk 2 y(u) = (5.14) 
The quantity of interest is the Wronskian between the two solutions defined by 



it — > oo 



y x {u,k) = ^ + 0(u- 
y 2 {u,k) = 1 + 0(u - 1) it — > 1 



(5.15) 



To perform a WKB analysis we change variables x = y/u and \ = V x 7 — %U- The 
resulting equation is of the Schrodinger type 



- X "(x) + V(x) X = 



V(x) 



1 35x u - 70x b 



ix 3 k 2 



4 x 2 (x 6 - l) 2 

The WKB analysis is completely similar to the previous case. The result is 

w(k) 



(5.16) 
(5.17) 



W( yi (u),y 2 (u)) 
w{k) 



u A — u 
6^ 



sin(a/c) 



with 



a 



dy 



%k 3 

v^£(|) 
" 3 r(|) 



(5.18) 
(5.19) 

(5.20) 



There is also a pair of extra poles at k = ±ir/a. This can be cured replacing the 
Wronskian by 

6^3 



w(k) - 
K J nk k 2 

From where we read the glueball masses as 



sin(afc) 



m r 



—n, n = 2,3, 
a 



(5.21) 



(5.22) 



in good agreement with the results of 0]. Analytically continuing to k — > ik the 
desired Wronskian is 

w(k) = ^-^ T sinh(ak) (5.23) 
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